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ABSTRACT 

To understand magnetic diffusion, momentum transport, and mixing in the 
interior of the sun, we consider an ideahzed model of the tachocline, namely 
magnetohydrodynamics (MHD) turbulence on a /3 plane subject to a large scale 
shear (provided by the latitudinal differential rotation). This model enables us 
to self-consistently derive the influence of shear, Rossby and Alfven waves on the 
transport properties of turbulence. In the strong magnetic field regime, we find 
that the turbulent viscosity and diffusivity are reduced by magnetic fields only, 
similarly to the two-dimensional MHD case (without Rossby waves) . In the weak 
magnetic field regime, we find a crossover scale {Ljij from a Alfven dominated 
regime (on small scales) to a Rossby dominated regime (on large scales). For 
parameter values typical of the tachocline, is larger that the solar radius so 
that Rossby waves are unlikely to play an important role in the transport of 
magnetic field and angular momentum. This is mainly due to the enhancement 
of magnetic back-reaction by shearing which efficiently generates small scales, 
thus strong currents. 

Subject headings: Turbulence - Sun: interior - Sun: rotation 



1. Introduction 

Data from global helioseismology (Thompson et al. 2003) have shed some light on the 
internal rotation of the sun. Throughout the convective envelope, the rotation rate decreases 
monotonically toward the poles. Near the base of the convection zone, there is a sharp tran- 
sition between differential rotation in the convective envelope and nearly uniform rotation 
in the radiative interior. This transition region has become known as the solar tachocline 
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(Spiegel 1972), and its thickness has been constrained to be only a few percent of the solar 
radius by helioseismic data (Charbonneau et al. 1999). The rotation rate of the radiative 
interior is intermediate between the equatorial and polar regions of the convection zone. 
Thus, the radial angular velocity gradient across the tachocline is positive at low latitudes 
and negative at high latitudes. The tachocline is most probably located below the convec- 
tion zone (Charbonneau et al. 1999), but part of it might be included in the overshooting 
region where turbulent plumes coming from the convection zone can induce turbulent mo- 
tions. This overshooting layer contains a prominent latitudinal differential rotation. It is 
thus of primary importance to understand the physics pertinent to such a turbulent layer in 
presence of strong shear. In particular, as the tachocline links two regions of very different 
transport properties, theoretical predictions of transport coefficients are essential to under- 
stand its long term dynamics (e.g. its thickness). For instance, different mechanisms have 
been invoked to stabilize this transition layer against the radiative broadening during the 
solar evolution. First, since the radial mixing is ineffective due to the stable stratification 
in the tachocline, Spiegel & Zahn (1992) suggested that the efficient turbulent transport of 
angular momentum (eddy- viscosity) in the horizontal direction could be responsible for the 
confinement of the tachocline. Other models rely on the existence of a magnetic field fully 
contained in the radiative zone (Riidiger & Kichatinov 1996; Gough & Mclntyre 1998; Mac- 
Gregor & Charbonneau 1999). Specifically, Riidiger & Kichatinov (1996) have shown that 
this configuration could explain the dynamics of the tachocline if the effective viscosity is 
larger by at least 4 orders of magnitude than its molecular value. Gough & Mclntyre (1998) 
and MacGregor & Charbonneau (1999) have estimated the thickness of the tachocline by 
balancing the shearing of the poloidal field with the resistive dissipation of toroidal magnetic 
field. However, in this estimate, the turbulent diffusivity and viscosity might be more relevant 
than the molecular values. The tachocline is also important for the generation of magnetic 
field inside the sun (by means of a dynamo process), whose operation crucially depends on 
the values of effective (eddy) viscosity and magnetic diffusivity (e.g. Parker 1993). 

The purpose of this paper is to provide a consistent theory of momentum transport and 
magnetic field diffusion in the tachocline. In view of the tendency of two-dimensional (2D) 
turbulence in a stably stratified layer, we focus on the dynamics in the (local) horizontal 
plane orthogonal to density gradient, by incorporating the effects of rotation and toroidal 
magnetic fields. Specifically, the latitudinal differential rotation is represented by a large- 
scale shear fiow while the latitudinal variation in the Coriolis force is captured by using a 
local f3 plane. This model permits us to study the influence on turbulent transport of shear 
flow, Rossby and Alfven waves. The shear flow efficiently generates small scales by shearing, 
thereby enhancing the back-reaction of fluctuating magnetic fields. We show that in the 
limit of strong magnetic field, the turbulent viscosity and magnetic diffusivity are reduced 
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by magnetic fields only, with the same results as in the 2D magnetohydrodynamics (MHD) 
case (Kim & Dubrulle 2001). In the weak magnetic field regime, we find a crossover from a 
Alfven dominated regime (on small scales) to a Rossby dominated regime (on large scales). 
Using parameter values typical of the tachocline, we identify these two regimes and show that 
Rossby waves are unlikely to play an important role in the tachoclinic turbulent transport. 
This work complements our previous study of 3D hydrodynamical (HD) turbulence in the 
tachochne (Kim 2005; Leprovost & Kim 2006), within which the effects of shear flows (due 
to radial and latitudinal differential rotations) and the average rotation were investigated. 

The remainder of the paper is organised as follows: we provide our model of the 
tachocline and the governing equations in section 2, the results for the hydrodynamical 
(HD) and magnetohydrodynamical (MHD) cases in section 3 and 4 respectively. In section 
5, we discuss the implication of our findings for the tachocline. 

2. Model 

To model the tachocline, we consider the motion in a plane {y,z) independent of the 
third coordinate x. In the tachocline case, the x, y and z directions stand respectively for the 
radial, azimuthal and latitudinal directions. For simplicity, we assume that the Coriolis force 
varies linearly with latitude: CI — (Qq + Pz)'k. This model has originally been introduced by 
Rossby (1939) to study a homogeneous sheet of fiuid on a sphere, the effect of the sphericity 
being only in the variation of the rotation rate with latitude. The term containing /3 is 
responsible for the appearance of the so-called Rossby waves. We also assume that, on large- 
scale, there are a linear shear fiow (due to latitudinal differential rotation), Uq = —{Az)y, 
and a uniform toroidal magnetic field, Bq = B^y, where ^ > and Bq > without loss 
of generality. This magnetic field supports Alfven waves in the system whereas a shear 
flow will favor the creation of small scale in the direction of the shear and thus enhance 
the dissipation (Kim & Dubrulle 2001; Kim & Diamond 2002; Kim 2006). The generation 
of small scales by a shear flow enhances the effect of Lorentz force via the formation of 
strong current. Consequently Alfven waves tend to become dynamically more important 
than Rossby waves, as shall be shown shortly. 

Owing to the geometry of the problem, it is convenient to work with the vorticity of 
the fluctuations a; = (V x v) • x and their magnetic potential deflned by b = V x (ax). We 
study these fluctuations in the framework of the quasi-linear approximation (Moffatt 1978) 
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where the interaction between the small scale fields can be neglected: 

{dt + Uody)u = -pv,-{Body)V^a + uV^uj + U, (1) 
{dt + Uody)a = -VzBq + r/V^a . 

Here, we assume that the turbulence is driven mainly hydrodynamically (no force in the 
equation for the magnetic field) . For simplicity, we assume an unit magnetic Prandtl number 
{ly — rj) and introduce a time dependent Fourier transform to capture the shearing effect: 



where F stands for u, a or fy^ and the azimuthal wave number evolves as: dtk^ = kyA. By 
using the transformation: F = Fexp[i/(/c^t + kl/3kyA)] to absorb the diffusion term and by 
transforming the time variable from t to t — kz{t)/ky = kz{to)/ky + A{t — to), Eq. (1) can 
be rewritten: 



Here, a = j3 / Aky is the ratio between the rate of change in the rotation vector and the 
shearing rate and 7 = B^ky/ A is the ratio between the magnetic field and the shear. Note 
that the velocity and magnetic fields can be obtained from the vorticity and the magnetic 
potential by the identities: Vy — Vz = K = ^kzd and hz = —ikyd. Equation (2) 

can be solved with the initial conditions: cu — d — Oair — tq — kz{to)/ky. Note that the 
precise value of the initial conditions is irrelevant as we are mainly interested in the long 
time limit where the memory of the initial conditions has been lost. 

We are also interested in the transport of chemical species whose concentration is again 
split in a large-scale and a fiuctuating part: N — No + n. In the quasi-linear approximation, 
by assuming the molecular diffusivity of particle to be the same as the molecular viscosity 
(for the case of non-unit Prandtl number, see Kim 2006), we obtain the following equation 
for the ffuctuating density of particles: 

ky{l + r^) 



which is the same as the second equation in (2) up to a multiplicative constant. Thus the 
turbulent diffusivity of magnetic field is the same as the turbulent particle diffusivity. 
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3. Hydro dynamical results 

We first consider the case without a magnetic field (7 = 0) to elucidate the influence of 
Rossby waves on a sheared turbulence. Equation (2) with 7 = can be readily integrated 
with the following result for the fluctuating vorticity: 



'TO 



fT 

J ^igfctrctan r— arctan t] j 



To calculate the turbulence intensity and the turbulent transport of angular momentum 
{i.e. the eddy viscosity ut deflned as {vyVz) = —J^rdzUo — J^tA), we prescribe the statistical 
properties of the forcing fyj to be spatially homogeneous and temporally short correlated 
with the correlation time r^: 

(/«,(ki, h)Uik2, t2)) = Tf (27r)25(ki + ka) 5ih - ^2) ^(ka) . (5) 

Here, 0(k2) is the power spectrum of the forcing. We focus on the strong shear limit: 
^ = uky/ A -C 1, which is valid on characteristic scales of fluctuation L ^ Lg = 2n{h' / A)^^"^ . 
The calculation of the turbulence intensity and turbulent viscosity is similar to that of 
Kim (2005) which we refer the readers to for technical details. The following solutions are 
obtained: 



2(27r)2^2y ^2 • 

The results in Eq. (6) are the same as those in the 2D turbulence with /3 = (Kim & 
DubruUe 2001). That is, turbulence amphtude and transport are reduced only by shear 
stabilization, and Rossby waves have no influence on the intensity of turbulence or the 
transport of angular momentum. In the kinematic regime (7 = 0), we can also compute the 
magnetic diffusivity rjT, which enters the induction equation for the mean magnetic field and 
is defined as {aVz) — —rjTdzAo — —rjTBQ-. 



On 



^'^ ^ \i A2 I '^21 2 ^^^^ arctan - ~ ^ / d^k 0(k) sin^ a arctan - . (7) 



{2nyA^ J a^kl L a\ {2n)^p 



1 



1 



aJ 



Here a = P/Aky and a = kz/ky. Thus, in the kinematic regime, the magnetic field diffuses 
more slowly compared to the 2D HD case (the sine cardinal function being always smaller 
than 1). The last expression in Eq. (7) even shows that for large enough /3, the magnetic 
diffusivity is reduced by a factor Note that the same reduction apphes to the transport 
of particle as r]T — Dt- 
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4. MHD turbulence on a plane 

Combining the two equations of system (2), we can obtain the following equation for 
the magnetic potential: 

dr [(1 + r^)dra] - iadra + 7^(1 + T')a = -^/^ (8) 

The exact solution to Eq. (8) can be written in terms of Heun's confluent function; however, 
we here derive the physically more transparent results in the case of a weak and a strong 
magnetic fleld separately. The boundary between the weak (|7| ^ 1 where 7 = B^ky/A) 
and the strong magnetic regime (I7I -C 1) is given by = 2t{Bq/A. 



4.1. Wecik magnetic field 

For I7I ^ 1, we calculate the correction to the magnetic potential and vorticity up to 
second order in 7 as: 

Cj{t) = Cjq{t) + 7tDi(T) + 7^tD2(T) + . . . and aij) = + 7ai(T) + -f'^hir) + . . . , (9) 

where Ojq{t) is given by (4). Solving order by order, we find that Cj is an even function of 7 
whereas a is an odd function. The first correction to the hydrodynamical result is: 

ai(T) = ^^f^^ia[^rct^nr-^rct^nt] _ -^^ ^ (^q) 

Using Eqs. (10), we can easily calculate the second order (in 7) correction to the turbulence 
amplitude and turbulent stresses. Assuming that the forcing has the following symmetry 
property (f){ky) — (f){—ky), we find that the first correction to the turbulence amplitude 
vanishes as it is odd in ky. The turbulence amplitude is thus given by Eq. (6), showing that 
neither Alfven nor Rossby waves change the turbulence amplitude. On the contrary, the 
Reynolds and Maxwell stresses are modified as follows: 

, , Tf /" ,2i , /I X f 1 27^ . 2 111 

^"^"^^ ^ -(2^y^^^^^H2fc^ + ^^^M2"^^'""aJ/' 
/, , V ^^f f ,2. 7V(k) oia l]/3\2/3 

<w ~ -3(4:4 [2 ' 

Here, a — P/Aky, a — k^/ky and ^ = vky A. One can see that the correction to the mag- 
netic stress is much larger than that to the Reynolds stress (recall that ^ <^ 1). Comparing 
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the zeroth order solution for the Reynolds stress and second order for the Maxwell stress, 
we find that the former is dominant if 7^0;"^^"^/^ -C 1 (giving us the result of the 2D case). 
In the opposite case (7^0;"^,^^^/'^ ^ 1), the transport of angular momentum is mainly con- 
trolled by the magnetic part of the stress i. e. by Alfven wave turbulence. In particular, the 
turbulent viscosity Ut = {{vyVz) — {bybz))/A is now positive, reminiscent of the direct cascade 
of energy in MHD turbulence (Kim & DubruUe 2001). These results enable us to calculate 
the characteristic scale at which the transition from the Rossby wave turbulence (with neg- 
ative viscosity) to the Alfven turbulence (positive viscosity) occurs, Lji — 27r {AB^/P^uy/"^, 
which corresponds to the "Rhines scale" (referring to the scale which separates, in ordi- 
nary geostrophic turbulence, eddy turbulence from wave turbulence). Our result shows that 
Ln ~ rI^* ~ i?m^ as we assume a unit Prandtl number. However, it is easy to show that for 
f] ^ u (as is relevant for the tachocline) or equivalently for P„i ^ 1, the first correction to 
the magnetic potential is the same as previously. Consequently, the magnetic stress remains 
the same, with u being replaced by rj. Thus, for low Prandtl number, the crossover scale is 
Lu = 27i{ABq/ P'^riY^'^ ~ rU'^- This scaling is the same as that of Diamond et al. (2006) 
obtained in a slightly different context (MHD turbulence on a /3 plane without a shear flow). 

In a similar way, we can calculate the turbulent diffusivity 77^. To calculate the first 
two leading order terms, we need to keep terms up to the third order in 7 for the magnetic 
potential (not shown here for brevity). Here again, only the leading order term survives 
when assuming (f){ky) = (j){—ky) and thus the turbulent diffusivity remains the same as the 
one obtained in the hydrodynamical case [see Eq. (7)]. 



4.2. Strong magnetic field 

When I7I 3> 1, we seek for a WKB solution of Eq. (8). Up to second order in e = I/7, 
a solution for the homogeneous part can be found as: 



1 



1 



F(r)^A ^ exp i-arctanr-Fe^ij"(T) cos -G(t) . (11) 
^1_I_7-2L2 JLe J 

Here, A and (p are two integration constants: 

«W = ^ + t(t ^ + TT75)' ('2) 
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Using the method of variation of constants, one obtains the solution for the magnetic po- 
tential and the vorticity: 



a(r) = 



dt 



ifw{t) 



TO 



1 



exp 



arctan r — arctan t 



+e\H{T) - H{t)) 



sm-jGO 



la / 

TV 



(13) 



a;(r) 



To 



it ^'W^axp 
AG'{t)y/lT^ 

1 



la 



2 ^arctan r — arctan tj + e^(if(T) — i?(t)) 



X 



G\t) cos - {G{t) -G{t)] 



laG'jtf - 2r . 1. . 
^^^^-,^^esm-{G(r)-G(t)} 



We can then use Eqs. (13) to calculate the Reynolds and Maxwell stresses and the turbulent 
viscosity. In the strong shear hmit (^ ^ 1), we obtain the following: 



(14) 



Equation (14) shows that the turbulent viscosity is the same as in 2D MHD (Kim & Dubrulle 
2001) except for the additional term proportional to due to the transport reduction by 
Rossby waves. However, the strong magnetic field limit is only valid for ky <^ A/ Bq^ and 
thus the correction term is bounded from above: a <^ I3Bq/A?. Using typical solar values, 
So ~ 1 T, ^ ~ 6 X IQ-^ /3 ~ 2.7 x IQ-^^ m-^s~\ we obtain a < 0.75 and consequently 
the correction term (proportional to a^) is very small. 



The turbulent diffusivity can also be obtained in a similar way: 



2/3 



(15) 



which is the same as the 2D case. Therefore, both the turbulent viscosity and diffusivity are 
reduced by a strong large-scale magnetic field while the effect of Rossby waves is neghgible. 



5. Discussion 

In order to elucidate the turbulent transport in the tachocline, we considered a 2D 
model of turbulence on a /3 plane, in the presence of a large-scale latitudinal shear (with 
shearing rate .4) and a uniform toroidal magnetic field (intensity Bq). By using the quasi- 
linear approximation, we have computed the turbulent viscosity and diffusivity (of magnetic 
fields or particle) in the two limiting cases of weak and strong magnetic field. The intensity 
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of turbulence is not altered much by the presence of waves whereas the turbulent transport 
can be severely affected due to the effect of these waves on the phase of the field. 

In the case of a weak magnetic field, we have shown that there is a crossover scale 
Lji — 2t: [ABl / 13^ r]Y I from an Alfven dominated (on small scales) to a Rossby dominated 
turbulence (on scale large enough). In the Rossby dominated regime, we found a negative 
turbulent viscosity, suppressed only by the shear: vt oc —A~'^. In the Alvcn dominated 
regime, the turbulent viscosity is positive with the following scaling: Vt oc -Bq(/?^)^^. We 
also have shown that the turbulent diffusivity is not very much affected by a weak magnetic 
field and that for all scales, the turbulent diffusivity is suppressed only by Rossby waves: 
riT oc (the same as in the HD case). In the case of a strong magnetic field, we found 
that the turbulent viscosity and diffusivity are positive and the same as in the 2D MHD case 
(except for an additional factor in the turbulent viscosity which is very small in the tachocline 
context) and consequently are suppressed only by magnetic fields: r]T (x ut oc Bq^. 

These findings have interesting impfications for the dynamics of the tachocline. Using 
typical solar values, ~ 1 T, ^ ~ 6 x 10"^ s-^ /3 ~ 2.7 x 10"^^ m-^s-^ , u ~ IQ-^m^s-^ 
and T] ~ lm^s~^, we can calculate the boundary between the weak and the strong magnetic 
regime. Lb — 2t:Bq/A ~ 10^ m and the scale at which the transition from Alfven to Rossby 
turbulence occurs, Lr = 27r{ABQ/ P^r]Y^^ ~ 10^° m. Let us recall that this scale was 
calculated in the weak magnetic field case {L > Lb) which is consistent in the tachoclinic 
case as L^ > Lb- Figure 1 shows the relative positions of these two scales compared to the 
radius of the Sun. We also show on this figure the minimal scale Lg = 2tx{v / A)^^"^ ~ 10^ m 
for which our strong-shear approximation holds. As relevant motion in the tachocline are 

L,! - 10*^ 111 



Lr ~ 10^ m 



Lr ~ 10^° m 
i?o ~ 7 X 10^ m 



Fig. 1. — Sketch of the different relevant scales evidenced in this paper for tachochnic values 
of the parameters. Lr is the scale which separates the strong magnetic field regime from the 
weak one and Lb (obtained in the weak field approximation) is the scale which separates a 
Alfven dominated regime from a Rossby dominated one. We also included the smallest scale 
Ls for which the large shear limit is valid. 

probably on scales smaller than the solar radius, the dynamics of the tachocline is likely 
to be insensitive to the effect of Rossby waves. In the worst case, the dynamics would 
be that of the weak magnetic field with Alfven dominated turbulence (meaning a positive 
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turbulent viscosity). These results are mainly due to the enhancement of the back- reaction 
of fluctuating magnetic fields by shearing. 

The reduction of the turbulent diffusion and transport induced by Alfven waves may 
be problematic in the solar context, invalidating some of the mechanisms believed to be 
important for the tachocline dynamics (see also Diamond et al. 2006), e.g. the turbulent 
horizontal viscosity model of Spiegel & Zahn (1992). However, the important point in this 
particular scenario is that the horizontal viscosity can be much larger than the vertical one. 
We have shown (Leprovost & Kim 2006) that it is the case for a turbulence sheared more 
strongly in the radial direction than in the latitudinal one (as is the case in the tachocline 
where the radial shear is at least one order of magnitude greater than the latitudinal one). 
If the effect of the magnetic field is the same in the two directions (reduction by a factor 
Bq), this result should not be altered by the inclusion of a magnetic field. However, a detail 
analysis of this situation requires the study of a three-dimensional model and is outside the 
scope of this paper. 

Finally, we note that in this paper, we assumed the toroidal magnetic field to be uniform 
on a local Cartesian (/3) plane. However, on a sphere, there is a possibility of instability in 
presence of a band of toroidal magnetic field and latitudinal differential rotation (Oilman & 
Fox 1997). The turbulence arising from this instability can then be considered as a source 
of part of the forcing in our model. 

Wc thank D. W. Hughes and S. M. Tobias for useful comments. This work was supported 
by U.K. PPARC Grant No. PP/B501512/1. 
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